Instabilities of off-centered vortices in a Bose-Einstein condensate 
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We study numerically the excitations of off-centered vortices in a Bose-Einstein condensate. The 
displacement of a single vortex and the separation of a doubly quantized vortex are considered. We 
find that the core-localized excitations of the precessing vortices continue to feature the property 
which implies that the vortices are unstable. The core-localized, dipolar, and quadrupolar excitations 
are found to deform continuously as functions of vortex displacements and intervortex separation. 

PACS numbers: 03.75.Lm, 03.75.Kk, 67.40.Vs 
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I. INTRODUCTION 

Quantized vorticity which accompanies rotational mo- 
tion is a characteristic feature of superfluids. Vorticity in- 
volves singular lines, vortex lines which exist inside any 
closed path enclosing finite circulation. The change in 
the phase of the order-parameter field around the vortex 
hne is quantized in integer units of 27r. 

In a Bose-Einstein condensate (BEC) of an atomic 
gas , the first vortex was created by the JILA group . 
They used a two-component condensate of ^^Rb atoms 
and imprinted a phase winding of 27r onto one compo- 
nent. The position of the vortex core was clearly seen 
and also the precessional motion of the core was ob- 
served In a one-component condensate, the ENS 
group 1^ created a vortex state using a rotating trap 
which consists of an optical spoon and a magnetic trap. 
Condensates containing up to 4 vortices were observed. 
Their method can be understood in analogy with fluid 
motion in a rotating vessel. The idea of a rotating trap 
has been employed by many groups thereafter. 

In addition to these singly quantized vortices, multi- 
ply quantized vortices with phase windings 47r and Stt 
can also be formed Q , using a topological method intro- 
duced by Nakahara et al. [6j- However, the experimen- 
tally recorded image of the cloud of atoms does not show 
any indications for the splitting of multiply quantized 
vortices. This surprising apparent stability contradicts 
with what has been widely expected Therefore, it 

is urgent to understand the stability of these off-centered 
vortices and the multiply quantized vortices. 

The axisymmetric and rotationally symmetric (singly 
quantized) vortex states have been studied and a core- 
localized excitation with negative excitation energy is 
found H, 0, ^1 ■ Concerning this excitation, 

a relation to the instability of a vortex state [lol 
(through the transfer of the population from the con- 
densate mode to the core-localized excitation) and also 
to the precession motion 0, ^| of the vortex core have 
been pointed out. Nevertheless, the axisymmetric and 
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rotationally symmetric systems only mimic the experi- 
mental situations because (a) the axisymmetric vortex 
state can become a vortexfree state only through the off- 
centered vortex states, (b) a system with a vortex under 
precessional motion is no longer axisymmetric. 

The understanding of the excitation level with nega- 
tive energy involves a puzzling feature because the 
value is connected to both the direction and frequency 
of the precession motion and the instability of the vortex 
state. These two cannot be divided as long as we study 
the axisymmetric state. Therefore, investigations of the 
off-centered vortex state are necessary in order to under- 
stand both the instability and the precessional motion. 
The off-centered vortex state has been studied through 
the Gross-Pitaevskii equation [l7| and the Thomas- Fermi 
(TF) approximatio nfla. ITsf . But no solutions of the Bo- 
goliubov equations [l| have thus far been presented for a 
non-axisymmetric vortex state. 

The motivation for an analysis of non-axisymmetric 
doubly quantized vortex may be explained in a similar 
way: The instability of the doubly quantized vortex state 
has been pointed out for the rotationally symmetric vor- 
tex state It relates the existence of a mode with a 
complex excitation energy and its exponential increase 
of population. But the axisymmetric state only lasts 
for an infinitesimal time if the vortex is unstable. It is 
not known whether the system keeps having the feature 
(complex eigenvalue) on which the instability of doubly 
quantized vortex state depends. 

This paper presents a two-dimensional analysis of the 
excitation spectra supported by the condensate with an 
off-centered vortex, a doubly quantized vortex, and pairs 
of split vortices. We discuss the above questions of pre- 
cession, the instability of singly quantized vortex, and 
the doubly quantized vortex within a unified numerical 
framework. 



II. SINGLE OFF-CENTERED VORTEX 

We consider a two-dimensional {x, y) system in the ro- 
tating frame whose angular velocity is lo. The axis of 
rotation is perpendicular to the (cc, y) plane. Ideally, this 
represents not only the static condensate in the rotating 
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vessel but also the precessing vortex. The Hamiltonian 
is 

H{u:) = H -LJ-L (1) 
H = J ¥{-C\/^ + V)-^ + ^g¥¥Mdr (2) 

L = [i/\rx p)ifdr (3) 



where C = —h^/{2m), g = ATrh^a/m, and V is the confin- 
ing potential with a denoting the s-wave scattering length 
and m the mass. We use a; = (0,0, w) and r = {x,y,0). 
The time-dependent Gross-Pitaevskii (GP) equation may 
be written as 

{-CV^ + V{x,y)+g\^ix,y,t)\^ 

-u.rxp}4,{x,y,t) = ^n^^^^^,{4) 

where 4>{r) is the condensate wavefunction. The time- 
independent form is 

{-CV^ + V{x,y)-^l + g\^ix,y)\^ 

-u}-r X p}(j){x,y) = 0, (5) 

where fi is the chemical potential. We employ an 
axisymmetric harmonic trapping potential V{x,y) — 
muj^^{x'^ + y^) where uiu- — 2tt x 200. The linear num- 
ber density of particles N — J \(j){r)\'^dr is fixed to 1 x 
lO^^m"^ in this paper. The mass m — 38.17 x lO"^'' kg 
and the scattering length a = 2.75 nm for Na atoms are 
employed. 

The angular momentum of the condensate per particle 



/ (/)*(r X p)(/)(ir 
h J |0(r)|2dr 



(6) 



is extensively utilized in the following discussion and plot- 
ted in several figures. A condensate with a centered 
vortex has L ~ 1. Systems with an off-centered vor- 
tex have L between and 1. We also use the TF radius 
Rtf — 6.79 /im as the unit of length. 



A. Static Analysis 

The time-dependent Bogoliubov equations 
du{x,y,t) 



ih 



dt 

- u; • r X p)u{x,y,t) 
- g(j>'^v{x,y,t), (7) 
dv{x,y,t) 
dt 

= - (-CV^ + V- ^1 + 2.g|(/.p + w • r X p) v{x, y, t) 
+ gcb*^u{x,y,t), (8) 



yield the excitation spectra of the condensate which fol- 
lows the time-dependent GP equation Eq. Q). They re- 
duce to the Bogoliubov equations 



u) ■ r X p) u 

a2„ 



{-CV^ + V - + 2g\(l)\^ 

-gcf^v = eu, (9) 
- (-CV^ + V -n + 2g\(j)\^ +uj-rxp)v 

+g(j)*'^u = ew.(lO) 

when the system is static. 

The criteria for the applicability of the Bogoliubov 
Eqs. ijnji and (|10|) . instead of the time-dependent Bogoli- 
ubov Eqs. (O and ((SJ, depends on how fast the wave- 
function of the condensate changes. We measure the am- 
plitude of the transformation using 



v 



(|(-CV2 + V + - u; • r X p - ^)<?!)|) 
T% ■ max(|0|) 



where the chemical potential is fixed here 0, [20l |. 
The numerator is a maximum of the left-hand side of 
Eqs. Q and ©. Therefore, in the framework of the 
time-dependent GP equation 



V'{t) = max 



d(j){x,y,t) 



dt 



/max(|0(a:,y,t)|). (12) 



For variations of the displacement of the vortex core r, 
the V has local minima for an axisymmetric vortex state, 
a vortexfree state, and an off-centered vortex state with 
certain displacement r for a choice of w. 

We use off-centered vortex states which satisfy 



V < 0.0005 



(13) 



in the following analysis of the Bogoliubov equations. 
Equations H12|l and (|13l) imply that the wavefunction (p 
varies slowly enough in the time scale of the harmonic 
trap [V < 0.0005 < wtr/(27r)]. Instead of the right-hand 
side of Eq. (|13|) , for low values of the angular momentum 
(L < 0.2) the right-hand side is rather relaxed to 0.03. 

Each of the relaxation processes, which is numerically 
equivalent to following the time-dependent GP equation 
with imaginary time, begins with the TF density pro- 
file as the initial configuration on which the phase of an 
off-centered vortex with various displacements r is "im- 
printed" through multiplication. In the beginning of a 
numerical simulation, the angular velocity u; is chosen. 
Then the development with fixed lu is performed without 
restricting the displacement. The vortex moves freely 
during the development and ends as an axisymmetric 
vortex, a vortexfree state, or a vortex state which sat- 
isfies Eq. H13|) . The latter form a set of cj) and uj val- 
ues which fulfill Eq. (|13|) . Figure ^a) shows some of 
the density profiles of the condensate. The displace- 
ment of the vortex core r and the angular momentum L 
are calculated through Eq. (|SJ| from each of the resulting 
wavefunctions d>. Therefore, L and r are functions of uj 
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FIG. 1: (a) Density of the condensate along the x-axis. The 
angular momenta L are 0.1277, 0.6854, and 0.9694, respec- 
tively, (b) The solid line represents the displacement of a 
vortex core vs. angular momentum of the condensate per par- 
ticle. The displacement vanishes for L — 1. The dotted line is 
the analytical estimate {1 — {r / Rtf)'^}'^ by Guilleumas fT3 . 
These two plots overlap in a wide range of displacements, 
< r < 0.8. (c) Displacement of the vortex core r vs. the 
normalized additional energy AE' of the condensate in the 
presence of an off-centered vortex. The solid line represents 
the numerical result. The dotted line is obtained by fitting 
{1 - (r/i^TF)^}^/^ taken from Ref. il, Eq. (49). 



through (f). The relation between the angular momentum 
and the displacement is plotted in Fig.mb). 

Figure ^c) displays the normalized additional energy 
due to the existence of an off-centered vortex, defined 
through 



AE' = {E- Eo)/{E, - Eo) 



(14) 



where Eq, Ei, and E are the energies of a vortexfree 
condensate, a condensate in the presence of a centered 
vortex, and with an off-centered vortex, respectively. If 
the condensate is vortexfree, AE' = 0. Off-centered vor- 
tex states have < AE' < 1 and the centered vortex 
state has AE' = 1. 



B. Excitation Spectra 

Excitations from a BEG have been interpreted in the 
framework of the Bogoliubov equations in Ref. • They 
play a significant role in determining vortex stability. 
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FIG. 2: (a) Excitation spectra of a Bose-Einstein conden- 
sate with a centered vortex, an off-centered vortex, and with- 
out a vortex. The vertical scale measures the excitation en- 
ergy, Slab. The horizontal axis indicates the angular mo- 
mentum qg. The crosses, triangles, squares, and the bul- 
lets correspond to the angular momenta of the condensate of 
L — 0,0.38,0.61, and 1, respectively, (b) Excitation energies 
of the dipole modes (labelled -1-1 and —1) and the quadrupole 
modes (-1-2, —2). (c) The angular momenta of the compo- 
nents u (solid line) and v (dotted line) of wavefunction for 
the dipole (-1-1, —1) and quadrupole modes (-1-2, —2). (d) 
Energy e of the lowest core excitation in the rotating frame 
is negative and it satisfies |e| <C 1. This is in accordance with 
the weak instability of the GP equation, Eq. 11311 . 



Nevertheless, these modes in off-centered vortices have 
not been analyzed so far. An excitation whose energy is e 
is expressed with two components, u and v, of wavefunc- 
tion each with different momenta. They are determined 
as the solutions to the eigensystem of Bogoliubov Eqs. © 



and (|10|l . The modes with negative 



')dr are 



ignored. Two of the lowest excitations are the condensate 
mode (with £ = 0,u = 4>,v = 4>*) and the core-localized 
mode. We ignore the mode with the smallest \e\ as the 
condensate mode. 

The wavefunction expressed with u and v are obtained 
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through the Bogohubov equations Q and (|10|l . We de- 
note their angular momenta using 



A{u) 
A{v) 



_ J u*{r X p)udr 
fi J juprfr ' 
_ J v*{r X p)v dr 



(15) 
(16) 



These quantities correspond to the angular momentum 
of the condensate L defined in Eq. JHJl. Some algebra 
with the help of Eqs. (0, ©, and shows that the 
excitation energies e depend linearly on uj. We presume 
that the wavefunctions 4>, m, and v do not depend on uj 
and energies fi and e do not have imaginary part. The 
coefficient of proportionality is 

_ Re [{A{u) -L}J \u\^dr + {Aiv) + L} J \v\^dr] 
" /(|u|2 + |w|2)rfr ■ 

(17) 

We introduce an excitation energy in the laboratory 
frame 



(18) 
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Equation (|17|l includes A{u), A{v), and L because an ex- 
citation level has two components u and v of wavefunc- 
tion and their phases are always affected by the phase of 
the condensate. For axisymmetric systems, the qa def- 
inition reduces to an integer quantum number llClj |. 
The excitation energies e and eiab are normalized by the 
trap unit fiuott from here on. Figure \^&) displays the 
computed excitation energies, e, and the corresponding 
angular momenta, qg. Some of the modes, for exam- 
ple those with (eiabi^e) = (2,0),(1,±1) maintain simi- 
lar values of qg and eiab upon variations of L. We call 
the modes at (2, 0) and (1, ±1) the breathing H and the 
dipole modes, respectively. Some of the core-localized 
modes have (eiabj^e) = (0,0). These are discussed in 
Sec, urn The modes with {eub,qe) ^ (1, ±1), (1.5, ±2) 
are discussed in Sec. IIIDI 



C. Core-Localized Mode 

For any value of L, the system has one core-localized 
mode. The modes always displays a sharp peak in the 
amplitude |up + |wp of the wavefunction at the vortex 
core. The core-localized modes still feature a major dif- 
ference between the system with L = 1 and those with 
< L < 1. 

When L = 1, the vortex is in the center of the system. 
The angular momentum qg — —1 and the excitation en- 
ergy eiab is negative [the lowest bullet ai qg — —1 in 
Fig. Ufa)]. It has a sharp peak in the vortex core. The 
sharp localization, the angular momentum, and the neg- 
ative energy are equivalent with those for axisymmetric 
calculations [9lll0l|. 



FIG. 3: Wavefunctions \uf of (a) the core-localized mode, (b) 
dipole mode with negative qe, (c) dipole mode with positive 
qe , (d) quadrupole mode with negative qe , and (e) quadrupole 
mode with positive qe- Solid lines denote contours of density. 
The core-localized mode has a peak of the amplitude in the 
vortex core. The peak in the quadrupole mode with negative 
qe is not seen in axisymmetric systems. The angular momen- 
tum L of the condensate is 0.531. 

Once the vortex becomes off-centered (0 < L < 1), 
both of the components u and v of the wavefunction lo- 
calize in the vortex core [Fig. Ola)]. Each of them fea- 
tures a sharp peak in the vortex core and its skirt extends 
towards the surface of the condensate. The angular mo- 
menta A{u) of the core-localized modes vary between 1 
and 4. The corresponding wavefunctions v have oppo- 
site angular momenta. Their average qg almost vanishes 
{\qg\ <^ 1); the points at (0,0) in Fig. |2Ia) mean this. 
The excitation energy plotted in Fig. El^d) remains below 
zero (—0.01 < e < 0) even in the rotating frame. 

The core-localized mode behaves differently between 
the system with L = 1 and that with < L < 1. One 
possible cause is the static approximation of the time- 
dependent Bogoliubov equations introduced in Sec. Ill Al 
While the system with L = 1 is treated with the Bogoli- 
ubov equations, systems with < i < 1 are described 
within the time-dependent Bogoliubov equations. 

D. Dipole and Quadrupole Modes 

Excitations having the lowest positive energy at qg = 
±1 and ~ ±2 are classified as the dipole and quadrupole 
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modes. Figure ^h) shows the energy for these modes. 
The energies of the dipole modes almost equal one trap 
unit (htoti)- The dipole modes are known to be respon- 
sible for the center-of-mass motion '2^. Therefore, they 
are not affected by various profiles inside the conden- 
sate. The angular momenta A{u) and A{v) are plotted 
in Fig. [21[c). These excitations show continual transfor- 
mation as a function of variation of the vortex displace- 
ment. 

The condensate in a two-dimensional harmonic poten- 
tial has two quadrupole modes. The splitting between the 
quadrupole frequencies eiab(+2) — eiab(— 2) is inverted at 
L < 0.4, while it is proportional to L in systems with 
a centered vortex |23|. The corresponding eigenoscilla- 
tion of the condensate, called the scissors mode, has been 
observed experimentally [24] to measure an angular mo- 
mentum of the condensate. The inversion may affect the 
observations of the angular momentum. 

Guilleumas et al. |l8l | estimate the splitting using sev- 
eral methods within the TF limit. Their Eqs. (33) and 
(42) and our numerical results in Fig. 2(b) agree well 
when the vortex is close to the axis (r/i?TF < 0.2). Their 
estimates stay positive and approach zero as the displace- 
ment of the vortex approaches i?TF ■ But the sign is in- 
verted in our numerical results. This is caused by an 
increase of the energy eiab(— 2) at lower L, depicted in 
Fig.lHb). The wavefunction u of the —2 mode displays a 
sharp peak at the vortex center; for example, in Fig.lJl^d). 
This peak is significant for lower L and not seen in the 
axisymmetric (non)vortex states. We presume that the 
existence of this peak is pushing up the energy eiab(— 2). 

For u! > 0.35wtr (corresponding to L < 0.2), the exci- 
tation energies of the modes with qg > 2 approach each 
other, and their wavefunctions become numerically indis- 
tinguishable from one another. The angular momenta of 
the excitations in this region {L < 0.2) are omitted in 
Figs. El for the sake of clarity. 



III. TWO VORTICES 

A vortex with the winding number two has been 
formed in a condensate of Na atoms Q . Multiply quan- 
tized vortices are energetically unfavorable against the 
array of the singly quantized vortices (e.g. Sec. 9.2.1 in 
Ref. 0), while the experimental data does not indicate 
any signs of the splitting. In this section, we investigate 
the excitation spectra of systems with two vortices to 
inspect their stability in the splitting process. 

Each of the GP calculations starts from an initial wave- 
function with two slightly displaced vortices under a 
fixed angular velocity w, which acts as an external ad- 
justable parameter. The initial wavefunction is formed 
from the TF density profile in the absence of a vortex, 
but multiplied with the phase factor e** around each of 
the initial vortex positions. The resulting angular mo- 
menta are controlled through this lu. The GP equation 
has a static solution for about 60 < w/(27r) < 110. 
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FIG. 4: (a) Density of condensate along x-axis. A condensate 
has two vortices and both of them are on the x-axis. Angu- 
lar momenta L are 1.198, 1.596, and 1.791, respectively, (b) 
Displacement of vortex core vs. angular momentum of the 
condensate par particle. The displacement will be zero when 
the angular momentum of the condensate is 2. The dotted 
line is the analytical estimate 2{1 — {r / Rtf)^}'^ taken from 
Eq. (24) in Ref. [3- The two plots agree well. 



The corresponding range of angular momenta is from 
1.2 to 1.8. Many vortices enter into the condensate 
for higher Loji^-K) > 110 and the vortices disappear for 
L>:l{2ix) < 60. Figure 0] shows density profiles of the con- 
densate for various angular momenta. 

Using the condensate wavefunctions, the excitations 
from it are calculated with the Bogoliubov equations, 
Eqs. © and ijini). We found that one or two conjugate 
pairs of modes with complex eigenvalues always exists in 
the calculated range of L. 

Figure 13d) shows the energy levels for two excitations 
with the lowest real part in the rotating frame. For 1.2 < 
L < 1.57, there occur two conjugate pairs of complex 
eigenvalues, both of them have a large imaginary part. 
For larger L, only one of them has a large imaginary 
part. Figure displays the wavefunctions of these lowest 
modes. One of them has an imaginary eigenvalue, while 
the other one has a real eigenvalue. Both of them display 
a peak in the vortex core. 

A doubly quantized vortex in a rotationally symmetric 
system has two negative excitations, those for qg = —1 
and —2. One of them with qg = —2 sometimes has a 
complex energy, depending on the particle density in the 
system, see Ref. 7J. This feature is equivalent to those 
of our results depicted in Fig. [S{d) for higher L. There- 
fore, the splitted vortices keep on having unstable nature 
which exist before the splitting 01 . 

Figure [SJa) shows that the modes with positive en- 
ergy stay with similar values of eiab and qg as functions 
of the separation of the vortices. The lowest positive- 
energy modes, with qg ~ ±1 and qg ~ ±2, are classified 
as the dipole and quadrupole modes. Their computed 
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FIG. 5: (a) Excitation spectra of a Bose-Einstein condensate 
with a pair of off-centered vortices. The vertical axis denotes 
the excitation energy, eiab. The modes with complex ener- 
gies are not plotted here because the definitions of eiab and 
qg in Eqs. 1171 and 11811 do not support modes with complex 
energies. The horizontal axis is the angular momentum, qg. 
The crosses, triangles, squares, and bullets correspond to the 
angular momenta of the condensate 1.198, 1.449, 1.569, and 
1.791, respectively, (b) Energy levels eiab of the quadrupole 
and dipole modes in the laboratory frame, (c) Angular mo- 
menta of the quadrupole and dipole modes. Solid lines show 
those of u, while the dashed lines correspond to v. (d) Two 
excitations with the lowest Re(e). Conjugate modes with neg- 
ative imaginary parts are not plotted. A system has two con- 
jugate pairs of complex eigenvalues or one pair of complex 
eigenvalues and one real eigenvalue. The real and imaginary 
parts are plotted independently. The points at L = 2 in (b) 
and (c) are taken from axisymmetric calculations. 




FIG. 6: (a) Wavefunction \u\ of a complex mode; the wave- 
function V of the complex mode is the conjugate of u. The 
energy levels of these modes are ±0.0241i. (b) (c) Wavefunc- 
tions |wp and \v\'^ of the lowest mode. The energy eigenvalue 



of this mode is eiab = - 
the condensate is 1.750. 



-0.2329. The angular momentum of 



IV. PRECESSION FREQUENCIES 



The above results are calculated in the rotating frame. 
It is not self-evident whether they may be interpreted as 
results for precessing vortices in the laboratory frame. To 
estimate the validity of the condensate wavefunction, we 
carry out time-dependent Gross-Pitaevskii calculations 
in the laboratory frame. The above results of the static 
GP equation are allowed to evolve in time according to 
Eq. 



angular momenta and energies in the laboratory frame 
are presented in Fig. Efb) and EJc). The modes near 
{qe,£iah) — (0,2) are breathing modes. The consistent 
behavior of the dipole, quadrupole, and breathing modes 
throughout the results in systems with an off-centered 
vortex [Figs.|2Ia)-|2Ic)], two vortices [Figs.|SIa)-[SIc)], and 
an axisymmetric vortex (L = 1, 2) proves the validity of 
our numerical procedure. 



The condensate shows precessional motion of the vor- 
tices for both the single- vortex and the two- vortex cases. 
Figure [7| shows the precessional angular velocity, Wpr. 
The rotation frequency uj that we obtained in the static 
calculations in Sees. Hll and IIIII are also plotted. The two 
frequencies agree well. It means that our results in the 
rotating frame can also be identified as those in the lab- 
oratory frame. The excitation modes, including the core 
mode, precess following the vortex core there. 
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FIG. 7: Circles represent the precessional angular velocities 
according to time-dependent calculations; lines show the an- 
gular velocity oj/tJtr that gives each value of angular momen- 
tum L. The vertical axis is normalized by the trap frequency, 
(a) Systems with one vortex, (b) Two vortices. 

V. DISCUSSION 

The excitation spectra of Bose Einstein condensates 
with one and two off-centered vortices has been com- 
puted. The excitation with negative energy and local- 
ization at the vortex core always exists not only in the 
axisymmetric state but also for the off-centered vortex 
configurations. We depicted the shape of wavefunctions 
and excitation energy of the core-localized, dipole and 
quadrupole modes in the off-centered vortex states with 
both one and two vortices. The invariance of dipole 
and breathing modes in the nonaxisymmctric cases shows 
that usage of of the Bogoliubov equation instead of the 
time-dependent Bogoliubov equation is valid. 

The rotation frequencies agree well with the precession 



frequencies of vortices in the laboratory frame. It enables 
us to consider the two problems, precession of the vortices 
and the vortex instability concerning the negative exci- 
tation, independently. The precessing vortices pursue to 
have a core-localized excitation with negative excitation 
energy in the framework of Bogoliubov equations. There- 
fore, the condensate with one precessing vortex maintains 
the instability as long as the vortex exists. In order to 
pose an answer to the problem between the precess- 
ing mode and the instability, this work is to be extended 
to treat finite-temperature systems 0, . 

Concerning the splitting of the doubly quantized vor- 
tex, the existence of two core-localized excitations is con- 
firmed throughout the calculated range of L. The vortex 
pair always has an instability 0, [2J| because systems 
with a vortex pair always support core excitations with 
complex excitation energies, c.f. Fig.|3Jd). 

Throughout the analyses, the system displays two 
quadrupole excitations. The excitation energies of the 
quadrupole modes Figs.EJb) andlSfb) are useful because 
these excitation is used to measure the angular momenta 
of the condensate experimentally. 
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